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Abstract 

The production of antiprotons in p + A reactions is calculated in a mi- 
croscopic transport approach employing hadronic and string degrees of free- 
dom (HSD). It is found that the abundancies of antiprotons as observed by 
the E910 Collaboration in p + A reactions at 12.3 GeV/c as well as 17.5 
GeV/c can approximately be described on the basis of primary proton-nucleon 
and secondary meson-baryon production channels for all targets. The trans- 
port calculations demonstrate that the antiproton rapidity distributions for 
heavy targets are sensitive to the p (or hadron) formation time in the nuclear 
medium. Within our analysis the data from the E910 Collaboration are rea- 
sonably described with a formation time of 0.4 — 0.8 fm/c in the hadron rest 
frame. 
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1 Introduction 



Since the first observation of antiproton production in proton-nucleus [0, ^ ^ and 
nucleus-nucleus collisions [|, ^ the production mechanism has been quite a mat- 
ter of debate. Especially in nucleus-nucleus collisions at subthreshold energies tra- 
ditional cascade calculations, that employ free A^A^ production and pN annihilation 
cross sections, essentially fail in describing the high cross sections seen from 1.5-2.1 
A-GeV H^, |]. Thus multiparticle nucleon interactions |1^ have been suggested as 
a possible solution to this problem. On the other hand it has been pointed out that 
the quasi-particle properties of the nucleons and antinucleons might be important 
for the p production process which become more significant with increasing nuclear 
density [|1T], |1^, 0, |T^, |15|, |T6[. Relativistic transport calculations have indicated 
that all the low energy data from proton-nucleus and nucleus-nucleus collisions are 
compatible with attractive p self energies (at normal nuclear matter density po) in 
the order of -100 to -150 MeV |l^, |l^ |T^. However, it had been stressed at that 
time that the high antiproton yield might also be attributed to mesonic production 
channels [jl9[ since pp annihilation leads to multi-pion final states with an average 
abundancy of 5 pions . 

With new data coming up on antibaryon production from proton-nucleus and 
nucleus-nucleus collisions at the AGS 



m 



1 and SPS \26, 27, 28, 29 



the p, A enhancement factors seen experimentally were no longer that dramatic as 
at SIS energies, however, traditional cascade calculations employing free production 
and annihilation cross sections again were not able to reproduce the measured abun- 
dancies and spectra 0, |3^, Q especially for S, S and Q, Q. Here additional 



collective mechanisms in the entrance channel have been suggested such as color 

In another language this 
56|, a precursor phenomenon for the 



rope formation |^ or hot plasma droplet formation ^5 
has been addressed also as string fusion 
formation of a quark-gluon plasma (QGP). 

The intimate connection of antibaryon abundancies with the possible observa- 
tion of a new state of the strongly interacting hadronic matter, i.e. the quark- 
gluon plasma, has been often discussed since the early suggestion in Ref. that 
especially the enhanced yield of strange antibaryons - approximately in chemical 
equilibrium with the other hadronic states - should be a reliable indicator for a 
new state of matter. In fact, the data on strange baryon and antibaryon production 
from the NA49 and WA97 Collaborations show an approximate chemical equilibrium 
|]39| , 0, ^ with an enhancement of the Q~ , yield in central Pb + Pb collisions 
(per participant) relative to pBe collisions at the same invariant energy per nucleon 
by a factor ~ 15. As pointed out in Ref. the data on multi-strange antibaryons 
at the SPS seem compatible with a canonical ensemble in chemical equilibrium. 

In Ref. 1^ Rapp and Shuryak have proposed multi-meson production channels 
for pp pairs to describe the antiproton abundancies in central Pb + Pb collisions at 
the SPS. Later on, Greiner and Leupold have applied the same concept for the 
A production by a couple of mesons including a or (for the s quark). In fact, 
their estimates have been supported by the microscopic multi-particle calculations 
in Ref. [HHl employing detailed balance relations for 2 hadron ^ n hadron transi- 
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tions and vice versa. Thus, in case of relativistic nucleus-nucleus collisions at AGS 
and SPS energies the low effective antibaryon absorption seen experimentally can 
naturally be explained by the backward reaction channels, i.e. multi-meson fusion 
channels. 

The latter mechanism, however, does not apply to proton-nucleus reactions since 
there is no longer an approximately thermal 'meson bath' feeding the baryon- 
antibaryon production channel as pointed out in Ref. Nevertheless, the p 



yield measured by the E910 Collaboration m p + A reactions at 12.3 GeV/c and 
17.5 GeV/c was found to be surprisingly high ^5[. An analysis within the Glauber 



model showed that either the annihilation cross section aann should be small com- 
pared to the cross section in free space or the p formation time Tp should be a 
couple of fm/c when employing the free annihilation cross section The authors 
of the latter work conclude that the 'naive' first-collision model with subsequent 
annihilation of antiprotons is incompatible with their data supporting the earlier 
finding from Ref. , where previous data from p + A and A + A reactions had 
been analysed within a similar model. 

In this work we will address p production in proton-nucleus collisions at 12.3 
GeV/c and 17.5 GeV/c within the covariant HSD[] transport approach |^ which 
is based on the concept of string formation and decay at invariant collision energies 
a/s > 2.6 GeV for baryon-baryon collisions and y/s > 2.3 GeV for meson-baryon 
collisions, while the low-energy reactions are described by known vacuum cross sec- 
tions or resonance formation and decay cross sections. A single sensitive parameter 
in this approach is the string duration or hadron formation time Tp ^ 0.7-0.8 fm/c 
which has been fitted to the proton rapidity distribution in nucleus-nucleus collisions 
at SPS energies |^ and kept fixed furtheron. For more details we refer the reader 
to the review [O. 



2 Qualitative considerations 

The physical picture within this transport approach for p production in p + A reac- 
tions is as follows: The impinging proton might create a NN pair in a first collision 
at position Xi (cf. Fig. 1) and the prehadronic quark- ant iquark configuration trav- 
els essentially in forward direction by a time 'jtf - where 7 denotes the Lorentz 
dilatation factor - until the hadronic p emerges at position X2- The hadronic state 
then propagates further in the same direction and reacts with the remaining nu- 
cleons elastically or inelastically. In the latter phase the dominant channel will be 
annihilation with a cross section 

„ 50[mb] 

(yannWs) = , (1) 

Vrel 



where the relativistic relative velocity is given by 



J - Mf - M|)2 - 4Mf M| 

2E1E2 ^ ' 



^Hadron-String-Dynamics 
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which is evaluated in the lab. frame. As shown in Fig. 3 of Ref. |^ the expres- 
sion (|l|) well describes the experimental data |50| on p annihilation for relative 
momenta up to 2-3 GeV/c. 

In rapidity space the produced antibaryons will be localized around midrapidity 
t/m] fast antiprotons hadronize later than slow p's due to the larger 7 factor and in 
view of (0) also experience a smaller annihilation cross section. Thus antiprotons 
from first-chance collisions in heavy nuclei should result in rapidity distributions 
that are highly asymmetric with respect to midrapidity and a larger suppression 
should be observed for y < ym than for y > ym, if the target is located at y = 0. In 
fact, a simple look at the experimental p rapidity distribution from Fig. 2 of Ref. 
| 25[| does not support the first-collision picture. 

On the other hand, in Ref. [|l6l it was found that a dominant fraction of an- 
tiprotons emerge from secondary pion-nucleon reactions in p + A reactions at lower 
energies. Such secondary reaction channels might also be a reason for the obser- 
vations of Ref. |]2S[ since the p production cross section in vrA^ reactions close to 
threshold is large compared to the production cross section in pp collisions (cf. Fig. 
6 of Ref. |TB|). The dymical effects of such secondary production channels are easy 
to map out within the string/hadron dynamical picture: In a first pN collision a 
'meson' {qq pair) M is produced at position xi which hadronizes after time •Jm'Tf at 
position X2- Then it propagates on average for a distance Xm ~ {o'totPo)~^ ~ 2 fm (in 
the lab. frame of reference) before reacting with another nucleon. In this reaction 
(at position X3) it may also produce a pre- iViV pair which hadronizes after time 7x7? 
in position X4. Then the antinucleon propagates in the same direction and may be 
rescattered, annihilated or emerge to free space. The antibaryons from secondary - 
or even ternary - processes have a much shorter way to escape the nucleus without 
being annihilated than those from first-chance collisions. Furthermore, their initial 
rapidity distribution will be shifted closer to target rapidity because the primary 
mesons are produced symmetrically around ym and only the energetic mesons will 
give a sizeable contribution to the p production provided that they hadronize in the 
medium. 

As mentioned above we here adopt the HSD transport approach that has e.g. 
been tested in Ref. [51] for pion production and proton stopping in p + A reactions 
at AGS energies using rp = 0.8 fm/c. Since the description of the related data 



is very satisfactory (cf. Fig. 8 in Ref. |5T[]) we continue with aspects of meson 



formation and secondary interactions, that are relevant for antiproton production 
and propagation in p + A reactions. 

In order to achieve a quantitative idea about the time scales for meson production 
and propagation we show in Fig. 2 for p + Au at 12.3 GeV/c the time dependence 
of the 'formed' meson number (solid line) in comparison to the number of string- 
like qq excitations (dashed line) that convert to mesons on a time scale of Tp in 
their individual rest frame. The reference frame for the transport calculation is the 
nucleon-nucleon cms. (at yiab ~ 1.63) with a Lorentz 7-factor of ~ 2.67 for the 
projectile proton and target nucleus, respectively, such that the dimensions of the 
target in space - as seen by the hadrons - are squeezed by the factor 7. The sharp 
increase of the dashed line at t= denotes the beginning of the p + Au reaction 
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with the excitation of strings that may lead to the formation of A^A^ pairs, too. 
The latter also appear as hadronic states on a similar time scale as the mesons - 
approximately 2 fm/c later - and may be annihilated or rescattered on the residual 
target nucleons. We recall that only the formed hadronic states are allowed to 
interact with the nucleons of the target again and to produce secondary mesons or 
baryon-antibaryon pairs. Thus antiprotons or mesons, that move fast with respect 
to the target at rest, may only form in the vacuum and not interact in the nucleus 
again. In fact, this is approximately the case for a -Be-target at this energy, while 
the dominant fraction of hadrons are formed inside the Au-imdeus and reinteract 
again (see below). For completeness we note, that the increase of the meson number 
for t > 7 fm/c is essentially due to the decays of mesonic resonances {p,uj, (p, K*) 
that increase the net meson number. 

The distribution in the invariant collision energy y/s is shown in Fig. 3 for p + Au 
at 12.3 GeV/c for baryon-baryon (BB) collisions (solid line) in comparison to the 
reactions of formed mesons with nucleons {mB, dashed line). The high energy peak 
of the 55-distribution (around ^/s ~ 5 GeV) is due to primary pN collisions; it is 
smeared out substantially due to Fermi motion. The low energy peak of the BB- 
distribution reflects multiple scattering of baryons in the target that, however, do 
not contribute to pp production due to the threshold of ~ 3.75 GeV (solid arrow). 
The i/i distribution for interactions of formed mesons with baryons (dashed line) 
extends to about 4.5 GeV, which is well above the pp production threshold of ~ 
2.8 GeV (dashed arrow). The distribution in the collision number above threshold 
is substantially larger for i?i?-reactions than for meson-nucleon collisions, however, 
the distributions have to be 'folded' with the respective pp production cross sections. 
We recall that in vrA^ reactions close to threshold this cross section is large compared 
to the production cross section in pp collisions (cf. Fig. 6 of Ref. |T6|) such that the 
mB channels might well compete with the BB channels. 



3 Comparison to data 

It should be noted that for the present analysis the perturbative production scheme 



of Ref. |T6[ has been employed due to the low multiplicities for p production (~ 2 x 
10~^ for p + Au at 12.3 GeV/c). All cross sections are the same as in the latter work; 
the only modification introduced is a new parametrisation of the individual rapidity 
(y) and transverse momentum (px) dependence for antiprotons in the elementary pN 
and vrA^ production channels, since the kinematical regime is substantially different 
from the latter study and the final states can no longer be described accurately by 
4- or 3-body phase-space, respectively. To this aim high statistics event generation 



has been performed within the LUND string model |Q - as implemented in the 
HSD approach - for p production and the resulting two-dimensional spectra have 
been fitted by a Gaussian in rapidity and exponential in p'^{y). As demonstrated in 
Refs. [|17| pl[] the LUND model works rather well for this purpose in comparison to 
the available data |SD|, however, uncertainties in the order of 50% for the elementary 



cross sections cannot be ruled so far since there are no data available for invariant 
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energies close to threshold. 

With the input cross sections defined we can step on with a quantification of 
the physics addressed above. In Fig. 4 we show the p rapidity spectra for p + Be 
at 12.3 GeV/c in comparison to the experimental data from [^. The dashed line 
{BBprim., upper part) gives the p spectrum from BB production channels prior 
to final state interactions (FSI). The contribution from secondary mB collisions is 
not visible on a linear scale for the Be target. When including p rescattering as 
well as annihilation the dash-dotted line [BB final) is obtained for BB production 
channels. It is clearly seen that p annihilation is more pronounced for low rapidities 
- with the target at y = - than for y > ym ~ 1.63. Thus for p + Be the pN 
production channel clearly dominates and the comparison with the experimental 
rapidity distribution (lower part of Fig. 4) essentially tests the adequacy of the 
elementary differential p production cross section in pN collisions which appears to 
be underestimated slightly (by about 30%) in the present approach using the input 
from Ref. |TB[ and assuming that the data point for y = 1.1 is solid. One might 



rescale the elementary cross section accordingly, however, the focus of our present 
work is not on fitting data. On the other hand, antiprotons from a feeddown of A and 
S° might also be contained in the data while they are excluded in the calculation^. 
Thus, provided that the reactions are described properly in the HSD transport 
approach, the cross sections for all targets should be low by about 30%, too, in case 
the p yield is dominated by the BB channels (see below). 

The numerical results for p + Cu and p + Au at 12.3 GeV/ c are shown in Figs. 5 
and 6 adopting the same assignment of the lines as in Fig. 4. It is clearly seen that for 
the Cw-target the contribution from mB collisions - as displayed by the short dashed 
line (mBprim.) - becomes visible and that the annihilation of p's from the secondary 
interactions is rather small as expected from the schematic Fig. 1. On the other 
hand, in case of the Au-target the antiprotons from A^A^ channels are annihilated 
to a large extent especially for y < y-m and the surviving p's in the experimental 
rapidity window 1.0 < y < 2.0 p5| practically stem from the secondary production 
channels. Nevertheless, the experimental data are underestimated on average again 
by about 30% as in case of the -Be-target which might be due to improper elementary 
cross sections and/or still missing production channels in the transport approach. 

The relative antiproton absorption according to Figs. 4-6 strongly depends on 
rapidity due to a larger formation time Tp and lower annihilation cross section with 
increasing rapidity y. In case of p + Be there is practically no p annihilation for y > 
2 since the fast antinucleons only form in the vacuum. For the p + Cu system this 
threshold is shifted to ?/ > 2.6 while for the p + Au system the annihilation becomes 
negligible only for y > 3. The sensitivity to the formation time tf is demonstrated 
in the lower part of the Figs. 5 and 6, where also the final antiproton rapidity 
spectra are displayed for tf = 0.4 fm/c (dashed line) and 1.6 fm/c (dotted line), 
respectively. In case of the light i?e-target no sensitivity to tf was found within 
the numerical statistics of the transport calculations since there is practically no 
contribution from the mB channels. As can be extracted from Figs. 5 and 6 a 

^In Ref. Eal such a feeddown contribution was estimated to be less than 5%. 
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shorter formation time Tp leads to an enhancement of the secondary production 
channel by mB collisions, whereas a larger formation time (~ 1.6 fm/c) essentially 
suppresses the contribution from secondary meson induced reactions in line with Fig. 
1. In fact, the experimental spectra for Cu and Au appear to be better described 
by '^F = 0.4 fm/c for the elementary cross sections employed. 

The average p multiplicity in the rapidity interval 1 < y < 2 at 12.3 GeV/c is 
compared in Fig. 7 as a function of target mass A with the data from pS) for different 
Tp. As mentioned before, for the i?e-target the calculations for Tp = 0.4, 0.8 and 1.6 
fm/c give approximately the same result, which in comparison to the experimental 
data is down by about 30% just scratching the lower level of the experimental error 
bars. For the Cw-target both calculations for Tp = 0.8 fm/c and Tp = 1.6 fm/c 
fall down too low. In case of the Au-target only the calculation for Tp = 1.6 fm/c 
is out of range in comparison to the data and also incompatible with the shape of 
the dN/dy spectra such that a formation time in the order of 1.6 fm/c should be 
excluded. 

We, furthermore, compare our calculations to the system p + Au at 17.5 GeV/c 
where the E910 Collaboration has taken data, too In this case (Fig. 8) the 

most important production channel is due to primary pN inelastic collisions, how- 
ever, the surviving antiprotons dominantly appear above midrapidity Um ~ 1-81 
due to the strong absorption. Within the rapidity window 1 < y < 2 of the data 
the contribution from secondary mB channels is almost compatible (thin solid line, 
mB final, compared to the dot-dashed line, BB final). The sum of the contribu- 
tions is reasonably well in line with the experimental data from |2^ (lower part) 
for Tp =0.8 fm/c (solid line). A slightly better description is achieved with Tp = 
0.4 fm/c for y ^ 1-1.5 in line with the data at 12.3 GeV/c, however, now the p 
rapidity distribution is slightly overestimated for y ~ 2. In view of the fact, that the 
elementary antiproton production cross section is underestimated by about 30% in 
our calculations, this provides a lower limit for the formation time Tp. 



4 Summary 

In this work we have analyzed antiproton production in p + A reactions at AGS 
energies within the HSD transport approach to examine the apparent low p an- 
nihilation seen experimentally |2^, Whereas the rather high antibaryon yield 



in nucleus-nucleus collisions at AGS and SPS energies can be attributed to multi- 
meson fusion channels employing detailed balance 0, the experimental data 
in proton-nucleus collisions at AGS energies appear to be compatible with primary 
proton-nucleon and secondary meson-baryon production channels as also found in 



Ref. for lower energies. 

Detailed calculations employing different hadronic formation times rule out long 
formation times of the order of 1.6 fm/c (in the hadron rest frame) since such times 
strongly suppress the meson reinteraction rate. The comparison with the data from 
the E910 Collaboration instead favors formation times of 0.4-0.8 fm/c when adopting 
the free antibaryon annihilation cross section. The latter cross section might be 
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somewhat reduced in the nuclear medium, but it is not hkely to be down by factors 
~ 5 as quoted in the analysis of Ref. p5[. We have pointed out, furthermore, that 
especially the high momentum/rapidity part of the p spectrum, which dominantly 
arises from primary pN collisions, shows a sensitivity to the p formation time since 
'fast' antiprotons may only form in free space and no longer get annihilated in the 
medium. Unfortunately, this rapidity regime is out of the E910 acceptance such 
that further experimental studies will be necessary to provide a final answer to the 
issue of antiproton formation times. 

The authors like to acknowledge valuable discussions with C. Greiner. 
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Figure 1: Illustration of primary pp production (in xi), formation (in X2) and prop- 
agation in the nucleus (upper part). The lower part shows the primary production 

of a pre- meson {qq pair) in Xi, its formation to a hadron in X2, the propagation of 
the hadron in the nucleus and rescattering with another nucleon in X3. In the sec- 
ondary interaction also a pie-pp can be created which hadronizes in X4 and further 
propagates in the nucleus or in the vacuum (see text). 
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- p+Au, 12.3 GeV/c 




t [fm/c] 

Figure 2: Time evolution of the string-like qq excitations (dashed line) and the 
formed meson number in p + Au reactions at 12.3 GeV/c within the HSD approach 
employing the 'default' formation time Tp = 0.8 fm/c. The time is given in the A^A^ 
center-of-mass system. 
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Figure 3: The differential number of baryon-baryon collisions (solid line) and meson- 
baryon collisions (dashed line) in p + Au reactions at 12.3 GeV/c within the HSD 
approach employing the 'default' formation time Tp = 0.8 fm/c. The solid arrow 
denotes the threshold for pp production in BB collisions and the dashed arrow the 
threshold for mB collisions, respectively. 
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Figure 4: Upper part: Antiproton rapidity distribution in p + Be reactions at 12.3 
GeV/c within the HSD approach (dashed hne, BBprim) without any final state 
interactions. The infiuence of rescattering and p annihilation is demonstrated in 
comparison to the dash-dotted line [BB final), which shows the final antiproton 
rapidity distribution. Note that there is no longer p annihilation for y > 2 for 
the formation time Tp = 0.8 fm/c adopted in the calculations. The contribution 
from mB channels to antiproton production in this system is not visible on a linear 
scale. Lower part: Comparison of the calculated p rapidity distribution (for Tp = 
0.8 fm/c) with the experimental data from Ref. |2^. Within the statistics achieved 
in the transport calculations there is no change in the p rapidity distribution when 
changing Tp in the interval [0.4, 1.6] fm/c. 
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Figure 5: Upper part: Antiproton rapidity distribution in p + Cu reactions at 12.3 
GeV/c within the HSD approach for baryon-baryon production channels (dashed 
hne, BBprim) without any final state interactions. The infiuence of rescattering and 
p annihilation is demonstrated in comparison to the dash-dotted line [BB final), 
which shows the final antiproton rapidity distribution from BB channels for a for- 
mation time tf = 0.8 fm/c. The short dashed line displays the contribution from 
secondary meson-baryon {mB) reactions without FSI while the thin solid line corre- 
sponds to the final p rapidity distribution from mB reactions including rescattering 
and annihilation. Lower part: Comparison of the calculated final p rapidity distribu- 
tion from BB and mB collisions (solid line for tf = 0.8 fm/c) with the experimental 
data from Ref. The dashed and dotted lines demonstrate the results for tf = 
0.4 fm/c and 1.6 fm/c, respectively. 
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Figure 6: Upper part: Antiproton rapidity distribution in p + Au reactions at 12.3 
GeV/c within the HSD approach for baryon-baryon production channels (dashed 
hne, BBprim) without any final state interactions. The infiuence of rescattering and 
p annihilation is demonstrated in comparison to the dash-dotted line [BB final), 
which shows the final antiproton rapidity distribution from BB channels for a for- 
mation time tf = 0.8 fm/c. The short dashed line displays the contribution from 
secondary meson-baryon {mB) reactions without FSI while the thin solid line corre- 
sponds to the final p rapidity distribution from mB reactions including rescattering 
and annihilation. Lower part: Comparison of the calculated final p rapidity distribu- 
tion from BB and mB collisions (solid line for tf = 0.8 fm/c) with the experimental 
data from Ref. The dashed and dotted lines demonstrate the results for tf = 
0.4 fm/c and 1.6 fm/c, respectively. 
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Figure 7: Antiproton rapidity distribution in p + A reactions at 12.3 GeV/c within 
the HSD approach for different formation times Tp = 0.4, 0.8 and 1.6 fm/c for Be, 
Cu and Au targets. The experimental data (full triangles) have been adopted from 
Ref. I2l. 
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Figure 8: Upper part: Antiproton rapidity distribution in p + Au reactions at 17.5 
GeV/c within the HSD approach for baryon-baryon production channels (dashed 
hne, BBprim) without any final state interactions. The influence of rescattering and 
p annihilation is demonstrated in comparison to the dash-dotted line {BB final), 
which shows the flnal antiproton rapidity distribution from BB channels for a for- 
mation time Tp = 0.8 fm/c. The short dashed line displays the contribution from 
secondary meson-baryon (mB) reactions without FSI while the thin solid line corre- 
sponds to the flnal p rapidity distribution from mB reactions including rescattering 
and annihilation. Lower part: Comparison of the calculated flnal p rapidity distri- 
bution from BB and mB collisions (solid line: tf = 0.8 fm/c; dashed line: Tp = 0.4 
fm/c) with the experimental data from Ref. p5 |. 
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